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Abstract 

Thermalization of quark matter and antiquark matter is studied with quark- 
quark-antiquark as well as quark-antiquark-antiquark elastic scatterings. Squared 
amplitudes of qqq — > qqq and qqq — > qqq at order a A s are derived in perturbative 
QCD. Solved by a new technique, solutions of transport equations with the squared 
amplitudes indicate that the scatterings qqq — ► qqq and qqq — ► qqq shorten the 
thermalization time of quark matter and antiquark matter. It is emphasized that 
three-parton and other multi-parton scatterings become important at the high par- 
ton number density achieved in RHIC Au-Au collisions. 
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1. Introduction 



Establishing a thermal state is critical in the formation of a quark-gluon plasma [1,2]. 
Initial high-energy heavy- ion collisions create deconfined matter which is not in thermal 
and chemical equilibrium. The problem of quark-gluon matter thermalization was studied 
early in parton dynamics in Refs. [3,4]. The 2-to-2 and 2-to-3 gluon scatterings produce 
a gluon-matter thermalization time greater than 1 fm/c [5-7]. In accounting for the 
elliptic flow of hadrons observed in Au-Au collisions at the Relativistic Heavy Ion Collider 
(RHIC), recent hydrodynamic calculations come to the conclusion that the thermalization 
of quark-gluon matter is finished within the time of 1 fm/c [8-11]. The thermalization 
process is fascinating and needs to be understood urgently. We have proposed 3-to-3 as 
well as 2-to-2 quark elastic scatterings to study the thermalization of quark matter [12]. 
But a thermalization time of about 1.8 fm/c is obtained from the mechanism. How rapid 
thermalization of quark matter comes about from quark scatterings is still a difficult and 
complicated task to explain. In this letter we include quark-quark-antiquark and quark- 
antiquark-antiquark elastic scatterings in an effort to understand better the thermalization 
of both quark and antiquark matter independent of gluon matter. 

We present four new aspects. The first is to recognize the importance of multi-parton 
scatterings at a high parton number density; the second is to calculate perturbatively the 
squared amplitudes of quark-quark-antiquark elastic scatterings which subsequently give 
the squared amplitudes of quark-antiquark-antiquark elastic scatterings; the third is to 
solve transport equations with a new technique; the fourth is to provide a new form for the 
formula of the interaction range of three-parton scattering. The squared amplitudes are 
derived by Fortran codes that implement the QCD Feynman rules. The new technique 
deals with a multi-parton scattering in terms of a sphere that moves with a particle. 
These four aspects form an integral part of our study of the thermalization of quark and 
antiquark matter. 

2. Probability of multi-parton scatterings 

The rapid thermalization is related to high gluon number density [13-15] at which 
triple-gluon elastic scatterings lead to a short thermalization time [16]. The importance of 
three-gluon scatterings is implied and we can thus speculate on possible or even sizeable 
contributions of other multi-gluon scatterings. We therefore give a Monte Carlo test on 
the occurrence of multi-parton scatterings from a parton distribution that is anisotropic in 
momentum space and inhomogeneous in coordinate space, as created in initial collisions 
of two Lorentz-contracted gold nuclei. Such a parton distribution has been generated by 
HIJING [17] for central Au-Au collisions at ^/saw = 200 GeV and cast into the form [18] 
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where the gold nucleus radius is i?4=6.4 fm, and k±, y, t, z and r are transverse momen- 
tum, rapidity, time, coordinate and radius in the transverse direction, respectively. One 
thousand and five hundred gluons are created from the distribution f(k±, y, r, z, 0.2 fm/c) 
by the rejection method within —0.3 fm < z < 0.3 fm in the longitudinal direction and 
r < Ra in the transverse direction. This volume of partons corresponds to a parton 
number density of 19.4 fm~ 3 . At such high density, we have to examine the occurrence 
of multi-parton scatterings. We assume a scattering of two partons when the distance of 
the two partons is less than a given interaction range. If n partons are within a sphere of 
which the center is at the center of mass of the n partons and of which the radius equals 
the given interaction range, a scattering of the n partons is taken to occur. When the 
scattering of n partons with a certain value of n is counted, other multi-parton scatter- 
ings are excluded. Therefore, the maxima of the numbers of 2-parton, 3-parton, 4-parton, 
5-parton, 6-parton and 7-parton scatterings at a given time are 750, 500, 375, 300, 250 
and 214, respectively, which give a sum of 2389. The numbers of n-parton scatterings 
with n — 2, • • •, 7 at t — 0.2 fm/c are denoted each by m n and plotted in Fig. 1. When 
the interaction range increases, more scatterings of not only two but also seven partons 
happen. All the scattering numbers saturate at an interaction range of 0.6 fm except for 
the three-parton scattering number which saturates at a smaller interaction range. If we 
consider the practical case when 2-parton scattering, 3-parton scattering, etc. may hap- 
pen at the same time, we need to estimate the probability of n-parton scattering. A guide 
for this probability is provided by R n = m„/2389, the ratio of the n-parton scattering 
number obtained above to the total number 2389. This ratio is drawn in Fig. 2. For 
an interaction range larger than 0.15 fm, the three-parton scattering becomes important 
in comparison to the two-parton scattering. At an interaction range of 0.62 fm, which 
acts as the radius of a unit volume of sphere, the 2-parton scattering has the occurrence 
probability of 30 %, the 3-parton scattering 20 %, the 4-parton scattering 14.6 %, and 
even the 7-parton scattering has 7.5 %. The scattering numbers and the ratios shown 
in Figs. 1 and 2 lead to the high parton number density of 19.4 fm~ 3 which implies a 
substantial occurrence of multi-parton scatterings. 

3. Quark-quark-antiquark elastic scatterings 

In Figs. 3 and 4 ten diagrams are plotted to show the elastic scattering of quark- 
quark-antiquark at order a^. The eight diagrams in Fig. 3 are related to two successive 
gluon exchanges and the two diagrams in Fig. 4 involve a triple gluon coupling. If the 
two initial or final quarks have the same flavor, the exchange of the two quarks generates 
other diagrams which are not displayed in Figs. 3 and 4. For uuu — > uuu, we must take 
into account 38 diagrams of which 28 diagrams contain quark-antiquark annihilation and 
creation, and then confront a cumbersome derivation of the squared amplitude for the 
scattering q(p 1 ) + q(p 2 ) +<?(— P3) — ► q{Pa) + q(P5) + q(-Pg), where pi = (E u pi) is the quark 
four- momentum. 

As an example, the derivation of the squared amplitude of the last diagram in Fig. 
3 is described here. Denote by qi, g and A the momentum, color and space-time index of 
the gluon created by the annihilation of the initial quark and antiquark, by q 3 , h and a 
of the other gluon, and by q 2 the momentum of the quark propagator, respectively. With 
the Feynman rules given in Ref. [19], the first Fortran code is designed to construct the 
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spin- and color-summed squared amplitude 
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where T g etc. are the color SU(3) group generators, g s is the quark-gluon coupling constant 
and gg = 47ra; s , respectively. The Dirac gamma matrices 7a and j a in Mf^_ are indicated 
by 7a' and 7 CT / in Mp^. The calculation of the traces is lengthy and cumbersome, but 
the second Fortran code is designed to accomplish the task. The Fortran code produces 
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colors 

which agrees with the result derived by hand. Including the average over the spin and 
color states of the two initial quarks and the initial antiquark, the spin- and color-summed 
squared amplitude is 
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which displays the shortest expression among all the individually squared amplitudes of 
the diagrams in Figs. 3 and 4, and the nine independent Lorentz-invariant momentum 
variables s 12 = (pi+p 2 ) 2 , s 2 3 = (P2+P3) 2 , s 31 = (p 3 +Pi) 2 , u 15 = (pi-p 5 ) 2 , u 16 = (pi-p 6 ) 2 , 
U24 = (P2 ~ Pi) 2 , U26 = {P2 ~ Pe) 2 , U34 = (P3 ~ Pi) 2 and u 35 = (p 3 - p 5 ) 2 . Interference 
of amplitudes of different diagrams totals 1004 nonzero terms. These interference terms 
are derived with the Fortran codes. The spin- and color-averaged squared amplitude of 
every diagram and the interference terms of different diagrams are used in the transport 
equation to study the effects of quark-quark-antiquark elastic scatterings. 

4. Transport equation 

Quark matter is assumed to have equal up-quark and down-quark distributions. The 
same assumption applies to antiquark matter. The mutual dependence of the evolution of 
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quark matter and antiquark matter derives from the two-body and three-body scatterings 
of quarks and antiquarks. With the assumption that the quark distribution in quark 
matter is the same as the antiquark distribution in antiquark matter, the variation of the 
up-quark distribution f\ is described by the transport equation 

g Q r d 3 p 2 d 3 p 3 d 3 p 4 ( , 44 , s 
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where the degeneracy factor gQ = 6 and the velocity of a massless up-quark v± — 1. The 
distribution function /i is a function of the position rj, the momentum p- x and the time t. 
The 2-to-2 parton scatterings and the 3-to-3 scatterings are represented by the first and 
the second terms on the right-hand side of the above equation, respectively. Equations for 
down-quark, up-antiquark and down- antiquark are written in a similar way. The squared 
amplitudes for the 2-to-2 scatterings are the spin- and color-averaged squared amplitudes 
of order a 2 obtained in Refs. [20,21]. 

The squared amplitudes | M uuu ^ uuu | 2 , | M uud ^ uud | 2 , | -M udu ^udu | 2 and | -M udd ^udd 
were obtained in the study of triple-quark elastic scatterings [12]. The calculation of 
I •Muuiwuuu 1 2 requires the individually squared amplitude of every diagram in Figs. 3 
and 4 and their interference terms, and the other diagrams arising from the exchange 
of the two quarks are also included. The squared amplitude | -M uud ^ uu d | 2 is based on 
diagrams F_, F + and and the diagrams generated from the three diagrams via the 
exchange(s) of the two initial quarks and/or the two final quarks; | -M u du^udu 1 2 is based 
on F_, F + and F* and the diagrams generated from all the diagrams in Figs. 3 and 4 via 
the exchange of the left quark line and the right quark line; | -M u dd^ u dd 1 2 is based on 
all the diagrams in Figs. 3 and 4 and the two diagrams generated from F_ and F + via 
the exchange of both the left quark line and the right quark line. Since the calculation 
of I 7W U uu^uflu 1 2 is based on the same set of diagrams as | -M UU u-> U uu | 2 , I -Muuu^uuu | 2 is 
obtained from the expression of | .M UU u^uuu 1 2 by the replacements Si 2 — > s 23 , s 23 s i2, 
S31 -> s 3 i, u 12 -> M32, u 13 -> u 3 i, u 2 i -> u 2 3, u 2 3 -> u 2 i, u 3 i -> ui 3 and u 32 -> u 12 . Sim- 
ilarly, I A^ ufl d^uud I 2 is obtained from | -M udd ^ ud d | 2 , I M ud ^udn | 2 from | .Muda^uda | 2 
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and | A^ u dd^udd 1 2 from | -M uu d^ uu d | 2 - These 3-to-3 quark-quark-antiquark scatterings 
bring more complexity to the evolution of quark matter than the 3-to-3 quark scatterings 
alone. 



5. Numerical results from the particle sphere technique 



The squared amplitudes | M uuu ^ uuu | 2 , | -A/f uud ^ U ud | 2 , | -Mudu^udu | 2 , | M u dd- 

, A4 UU u-+uuu I ; I -Muud^uud | , | A^udu^udu | ; | -^udd— »udd I ? I •^uflu^uufl | ; | -^uud-^uud I > 

I -Mudu^udu | 2 an d I -Mudd-+udd | 2 in Eq. (4) are calculated at a s = 0.5 and the Coulomb 
exchange divergence that is encountered is removed by the use of a screening mass for- 
mulated in Refs. [22-24]. Gluon propagators in Feynman gauge are used in the squared 
amplitudes. A screening mass in these propagators leads to some gauge dependence. The 
transport equation is solved until momentum isotropy is established while up and down 
quarks and antiquarks each with number 250 are generated from the anisotropic parton 
momentum distribution (1) inside the volume of —0.3 fm < z < 0.3 fm and r < 6.4 fm. 

When the interaction range is treated as being finite, a parton interacts only with 
surrounding partons. We define a sphere for a particle which is located at the center 
of the sphere. Every particle has a sphere and all the spheres have the same radius. 
The sphere is called particle sphere. The amount of particles in a sphere varies with the 
radius. The amount of particles also changes from one particle sphere to another. We 
denote the maximum amount of particles by m ps . The particle sphere radius is not the 
interaction range. When the sphere radius equals 1.6 fm, m ps = 75 is found for 1000 
partons generated from the anisotropic parton momentum distribution. We only allow 
particles inside the sphere to scatter. But in the volume defined above this restriction 
still produces almost the same number of n-parton scatterings as that when no sphere is 
defined. Therefore, the restriction can produce accurate result in n-parton scatterings. 
Since the amount of particles inside a sphere is an order of magnitude lower than the 
total number of partons inside the volume, the operation of searching for scatterings is 
reduced considerably. For example, the total number of operations for searching 5-parton 
scatterings at time t = 0.2 fm/c is 1000mp S 3.2 x 10 10 with the use of particle spheres, 
which is dramatically smaller than the number 1000 5 = 10 15 without the use of particle 
spheres. Particle spheres form a basic ingredient in our simulation of parton scatterings. 

When a parton traverses inhomogeneous or evolved parton matter, the radius-fixed 
sphere of the parton contains different amounts of partons at different time. After a 
multi-parton scattering is finished, a particle sphere is determined for each final parton. 
Before or after a scattering event a parton possesses a different particle sphere. However, 
the particle sphere radius is always fixed at 1.6 fm independent of time. While parton 
matter expands, the amount of partons inside a particle sphere is reduced. 

We define a scattering of two partons when the two partons have the closest distance 
less than the interaction range of \J (72^2/^. The 2-to-2 scattering cross section 02^2 is 
calculated with the spin- and color-averaged squared amplitudes in Refs. [20,21] which 
are regulated by the screening mass /id, 



M 



uuu^uuu 



+udd 



qq->qq 



'qq->qq 



167TS 2 9 



s + 2 



S[S 



2u 



,, , + - (s + 2u D ) In 



(5) 



6 



for the elastic scatterings of two quarks or antiquarks with the same flavor, 
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for the elastic scatterings of one quark and one antiquark with the same flavor. 

Defining a 3-to-3 scattering event depends on the parton positions and flavors. If a 
three-parton scattering occurs, the three partons must be in a sphere of which the center 
is at the center-of-mass of the three partons and of which the radius r hs is [16] 
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where m = 1 if | M 3 ^ 3 | 2 =| -M ud u^udu | 2 , | -A^ ud a_ ud d | 2 or | >f ufl a^uud | 2 , m = 2 
if | M 3 ^ 3 | 2 =| M uu ^ uun | 2 or | -M uu a_ uu d I 2 , m = 4 if | .M 3 ^3 | 2 =| -M uQfl ^ uflfl | 2 
or | A^ u aa^udd | 2 - Since | ^3^3 | 2 depends on the nine Lorentz-invariant momentum 
variables, we replace the phase-space integration by dui 4: du24du 34: dui 5 du25du 35 . However, 
this change is postponed to the next section. 

The use of a finite interaction range in determining 2-to-2 and 3-to-3 scatterings 
breaks the locality of the transport equation and thus Lorentz covariance. But the particle 
subdivision technique [25,26] based on the transformation / — > /' = if can restore the 
Lorentz covariance. At £ = 20, a solution of Eq. (4) is taken as the average of 20 runs of 
Fortran code starting from different sets of partons generated from the distribution (1) 
at t — 0.2 fm/c. The solution of the transport equation at the time of the order of 1.75 
fm/c is shown in Fig. 5 for various angles by the dotted, dashed and dot-dashed curves. 
The curves overlap and can thus be fitted to the Jiittner distribution, 

><*> = iHT^y < 9 > 

where the temperature of quark matter T = 0.27 GeV and fugacity A = 0.31. We get a 
thermalization time of 1.55 fm/c. 

6. Interaction range of three-parton scatterings 

The radius of a sphere defined in Eq. (8) is Lorentz invariant and depends on si 2 , 
s 23 and s 31 . The integration over p 6 gives 

Since | ^3^3 | 2 is a function of the nine Lorentz-invariant momentum variables, we 
replace d 3 p± and d 3 p^ with duudu2idu 3 4 : dui^du2bdu 3 ^ 1 

dundu2idu 3 4 1 = J 4 dp4 X dp 4y dp4 Z , (11) 
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du 15 du 25 du 35 = J 5 dp 5x dp 5y dp 5z , 
where u u = (p 1 - p 4 ) 2 , u 25 = (j>2 - P5) 2 and 

J 4 = 8EiE 2 E 3 I (v 4 - vi) • [(v 4 - v 2 ) x (v 4 - v 3 )] |, 

J 5 = 8E 1 E 2 E 3 I (v 5 - vi) • [(v 5 - v 2 ) x (v 5 - v 3 )] |, 



(12) 
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where V; is the velocity of the ith parton, and in the center-of-momentum frame the 
energies of three initial partons in the 3-to-3 scatterings are 
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with s s = si 2 s 23 + s 23 s 3 \ + S31S12 and the total energy of the three partons E\ + E 2 + E 3 = 
y/s = \fs12 + s 23 + S31 is obtained from the definition s = (pi + p 2 + p 3 ) 2 ■ Integrating 
over -u 25 to remove 8{E\ + E 2 + E 3 — £4 — E 5 — E 6 ), so we get 
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where in the center-of-momentum frame the energies of the three final partons are 
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and the absolute value of the derivative of E\ + E 2 + E 3 — E 4 — E 5 — E 6 with respect to 
u 25 is 
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which is evaluated at u 25 that is determined by the energy conservation relation E x + E 2 + 
E 3 - E 4 - E 5 - E e = 0. Inside Eqs. (21) and (22), 

d = (v lx - v 2x )(v ly - v 3y ) - (v ly - v 2y )(v lx - v 3x ), 

W = 2(1 - vi ■ v 2 )(l - v 2 ■ v 3 )(l - v 3 • vi), 

Wi = Vi z (Vj - V k ) 2 + Vj z (vj - V k ) • (V k - Vi) + V kz (v k - Vj) • (Vj - Vi), 

Wii = W 2 - W [(v ix - v kx ) 2 + (v iy - v ky ) 2 ], 

W m = Wii + W (v ix Vky - VjyV kx ) 2 , 

W ij0 = + W (v ix v ky - v iy v kx )(v kx v iy - v ky v ix ), 

where the subscripts i,j and k allow the three cases: i = 1, j = 2, k = 3; i = 2, j = 3, 
k = 1; i = 3, j = 1, k = 2; and 
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with n = 4, 5. The quantity u 25 determined by energy conservation gives 

u 2e = -s 12 - s 23 - u 2A - u 25 ( 23) 

that appears in | ^3^3 | 2 . Replacing tt 14 with w 16 in w 14 = — (s 12 + s 31 + tt 15 + u 16 ), we 
obtain 

^ = 8(2^ / ^^16^24^34^35 " ( 24 ) 

The variables M15, Wi6, w 2 4, M34 and u 35 have different minima: 



Mlmin = -V( S 31 + S s){sj 2 + S a )/(s^ 3 + S s ) = -(s 31 + S i2 ), (25) 

for -u 15 and -u 16 , 

«2min = ~\/ (s? 2 + S s )(s| 3 + S s )/(s§i + S s ) = -(Si 2 + S 23 ), (26) 

for -u 24 and 

«3min = ~V (S 23 + S S )(s§l + S *)/( S 12 + S s) = ~( S 23 + S 3 l), (27) 

for M34 and M35. In the center-of-momentum frame E 4 + E 5 + Eq = y/s and P4+P5+P6 = 0. 
While the three vectors p 4 , p§ and p^ are parallel, \ p 4 \, \ P5 | or | p§ | takes its maximum 
a/s/2. If, for example, p\ and p 5 point to opposite directions and pi, p 4 and p*6 point 
to the same direction, w 15 by its definition reaches — 2i? 1A /s, i.e. its minimum. Hence, 
when one of variables U15, u±q, u 24 , u 34 and W35 takes its minimum, the three vectors p 4 , 
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p 5 and Pq and one of pi, p 2 and p 3 are parallel. The minima set useful bounds to the 
Lorentz-invariant momentum variables by 



u u + u 15 + u 16 = u lmin 



(28) 



«34 + U35 + U 3e = M 3min 
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(29) 
(30) 



with u 36 = (p 3 — Pe) 2 . Eq. (28) means that u 15 and u 16 cannot take the minimum Ui min 
at the same time. If u 15 = Ui m { n , U14 = u±q = 0. Eq. (28) is a plane equation in three- 
dimensional space. The distances of the -u^-intercept, Mis-intercept and -u 16 -intercept of 
the plane to the origin are the same, — w lmin . Since u u < 0, < and u 16 < 0, the 
regions of the three variables are the interior of the triangle in the plane of which the 
three sides connect the Mu-intercept, Mis-intercept and -Ui6-intercept. The triangle is an 
equilateral triangle that has a side length of —\^2ui min . Similar discussion applies to Eqs. 
(29) and (30). The three planes indicated by Eqs. (28)- (30) are parallel. 

7. Summary 

We have derived squared amplitudes for qqq — > qqq and qqq — > qqq at the tree level. 
Feynman diagrams depicted in Figs. 3 and 4 show the 3-to-3 scatterings which lead to 
the interplay of quark matter and antiquark matter. These squared amplitudes form new 
contributions in transport equations for quarks and antiquarks. We have presented the 
particle sphere technique to solve the transport equations while the formula for the in- 
teraction range of three partons is obtained by an integration over the Lorentz-invariant 
momentum variables in the equilateral triangle regions. It is shown by the solutions of the 
transport equations that the thermalization time of quark matter and antiquark matter 
is of the order of about 1.55 fm/c. This reduced thermalization time is an effect of quark- 
quark- antiquark and quark-antiquark-antiquark elastic scatterings that suggests the im- 
portance of multi-parton scatterings at high parton number density. Another interesting 
effect of three-body elastic scattering is shown on heavy quark momentum degradation in 
quark-gluon plasma [27]. 

Acknowledgements 

This work was supported in part by the National Natural Science Foundation of 
China under Grant No. 10675079. X.-M. thanks CM. Ko for bringing his interesting 
work to the author's attention during quark matter 2006 conference at Shanghai. 



10 



References 



[I] STAR Collaboration, J. Adams, Nucl. Phys. A757 (2005) 102. 
[2]PHENIX Collaboration, K. Adcox, Nucl. Phys. A757 (2005) 184. 
[3]E. Shuryak, Phys. Rev. Lett. 68 (1992) 3270. 

[4]K. Geiger, Phys. Rev. D46 (1992) 4965; 

K. Geiger, Phys. Rev. D46 (1992) 4986. 
[5]G.R. Shin, B. Miiller, J. Phys. G29 (2003) 2485. 
[6]Z. Xu, C. Greiner, Phys. Rev. C71 (2005) 064901. 
[7]S.M.H. Wong, Phys. Rev. C54 (1996) 2588. 

G.C. Nayak, A. Dumitru, L. McLerran, W. Greiner, Nucl. Phys. A687 (2001) 457. 
[8]P.F. Kolb, P. Huovinen, U. Heinz, H. Heiselberg, Phys. Lett. B500 (2001) 232. 

P. Huovinen, Nucl. Phys. A715 (20 03) 299c. 

[9]D. Teaney, J. Lauret, E.V. Shuryak, [rmcl-th/0110037 

E.V. Shuryak, Nucl. Phys. A715 (2003) 289c. 
[10]T. Hirano, Phys. Rev. C65 (2001) 011901. 

K. Morita, S. Muroya, C. Nonaka, T. Hirano, Phys. Rev. C66 (2002) 054904. 

[II] K.J. Eskola, H. Niemi, PV. Ruuskanen, S.S. Rasanen, Phys. Lett. B566 (2003) 187; 
K.J. Eskola, H. Niemi, PV. Ruuskanen, S.S. Rasanen, Nucl. Phys. A715 (2003) 561c. 

[12]X.-M. Xu, R. Peng, H.J. Weber, Phys. Lett. B629 (2005) 68. 
[13]K.J. Eskola, K. Kajantie, K. Tuominen, Phys. Lett. B497 (2001) 39. 
[14]M. Gyulassy, P. Levai, I. Vitev, Nucl. Phys. B594 (2001) 371; 

M. Gyulassy, I. Vitev, X.-N. Wang, P. Huovinen, Phys. Lett. B526 (2002) 301. 
[15]F. Cooper, E. Mottola, G.C. Nayak, Phys. Lett. B555 (2003) 181. 
[16]X.-M. Xu, Y. Sun, A.-Q. Chen, L. Zheng, Nucl. Phys. A744 (2004) 347. 
[17]X.-N. Wang, M. Gyulassy, Phys. Rev. D44 (1991) 3501; 

M. Gyulassy, X.-N. Wang, Comput. Phys. Commun. 83 (1994) 307; 

X.-N. Wang, Phys. Rep. 280 (1997) 287. 
[18]P. Levai, B. Miiller, X.-N. Wang, Phys. Rev. C51 (1995) 3326. 
[19]R.D. Field, Applications of Perturbative QCD, Addison- Wesley, Redwood City, 1989. 
[20]R. Cutler, D. Sivers, Phys. Rev. D17 (1978) 196. 
[21]B.L. Combridge, J. Kripfganz, J. Ranft, Phys. Lett. B70 (1977) 234. 
[22]T.S. Biro, B. Miiller, X.-N. Wang, Phys. Lett. B283 (1992) 171. 
[23]K.J. Eskola, B. Miiller, X.-N. Wang, Phys. Lett. B374 (1996) 20. 
[24] S. A. Bass, B. Miiller, D.K. Srivastava, Phys. Lett. B551 (2003) 277. 
[25]B. Zhang, M. Gyulassy, Y. Pang, Phys. Rev. C58 (1998) 1175. 
[26]D. Molnar, M. Gy ulassy, Nucl. Ph ys. A697 (2002) 495. 
[27] W. Liu, CM. Ko, |nucl-th/0603004] 



11 



800 — — ~ — --r 




Interaction Range (fm) 



Figure 1: Curves from top to bottom show the numbers of n-parton scatterings with 
n — 2, • • •, 7, respectively. 



0.4 




0.25 0.5 0.75 

Interaction Range (fm) 



Figure 2: Curves from top to bottom show the ratio R n with n — 2, • • •, 7, respectively. 
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Figure 3: Two-gluon-exchange induced scatterings of two quarks and one antiquark. 
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Figure 5: Quark distribution functions versus momentum in different directions while 
quark matter arrives at thermal equilibrium. The dotted, dashed and dot-dashed curves 
correspond to the angles relative to one incoming beam direction 9 = 0°, 45°, 90°, respec- 
tively. The solid curve represents the thermal distribution function. 
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